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HIGHLIGHTS 


• The co-combustion kinetics of beet/bituminous coal was studied. 

• The co-combustion kinetics of switchgrass/bituminous coal was studied. 

• The dominant mechanisms of co-combustion at pre-peak and post-peak were determined. 
. There were kinetic compensation effects for selected samples within wide heating rates. 
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The kinetics of co-combustion of two herbaceous biomass species (beetroot and switchgrass) with 
bituminous coal was evaluated using a non-isothermal thermogravimetric method (TGA). The effects 
of heating rates, sample kinds and blend ratios on the kinetics of co-combustion were revealed based 
on two-stage scheme, which referred to the two stages before and after the maximum combustion rate. 
The kinetic parameters and combustion reaction mechanisms were tested by combining the isoconver- 
sional method (Ozawa-Flynn-Wall) and the Coats-Redfern method in order to find out the kinetics 
characteristics responsible for the combustion of the samples. There were significant difference between 
combustion kinetics of bituminous coal and biomass. The blending ratios and the heating rates had have 
certain effect on the mechanisms of switchgrass/coal blends at pre-peak. The dominant mechanisms 
associated with co-combustion kinetics for beetroot/coal at pre-peak and post-peak were described by 
the Avrami-Erofeev equations. The mechanisms of co-combustion for switchgrass/coal samples at 
pre-peak were described by the Avrami-Erofeev equation, Z-L-T equation or Anti-Jander equation; 
however, at post-peak, their mechanisms were described by the Avrami-Erofeev equation. The 
general kinetic compensation effect correlations were deduced for all samples within the heating rates 
of 10-90 “C/min. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Biomass is a renewable energy source that will provide 
environmental and economic benefits, and biomass residues are 
readily available in large quantities [1]. A wide variety of biomass 
materials, including herbaceous and woody materials, agricultural 
residues and energy crops are used [2], Switchgrass, a perennial 
warm-season bunchgrass native to North America, where it occurs 
naturally from Canada southwards into the United States and Mex¬ 
ico, has also been identified as having potential as an energy crop for 
Eastern Canada and for the US [3] . It is desirable as a potential 
biofuel in direct combustion systems because of its excellent burn 
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quality, ease of management, and high yield [4], Beetroot residue 
is a by-product from beet sugar production, and it is abundant espe¬ 
cially in the US, but they are currently used for lower-value products 
[5], Burning biomass alone cannot give very high energy output 
compared with coal combustion because of its high volatile matter 
and moisture content [6-8], Co-combustion of biomass and coal in 
existing power plants is a most promising option for application 
with renewable fuels [9,10], Muthuraman et al. detect that reactiv¬ 
ity of coal had been improved upon blending with wood [11], Kinet¬ 
ics of co-combustion of coal and pine sawdust blends have also been 
studied [12], the results reveal that the first step of biomass oxida¬ 
tion and coal combustion is effectively described by first order reac¬ 
tion mechanisms, and the second step of biomass combustion is 
more influenced by diffusion mechanisms. Kinetics of coal, palm 
shell and their blends are determined using the Doyle’s and 
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Nomenclature 

a degree of conversion 

A the pre-exponential or frequency factor (s -1 ) 

dm/dt burning rate (mg min -1 ) 

E a apparent activation energy (J mol -1 ) 

M mass fraction of sample 

R the universal gas constant (J/(mol K)) 

t temperature (°C) 

T temperature (K) 


Greek letter 

P heating rate (°C min -1 /°C s -1 ) 

Subscripts and superscripts 
0 initial 

oo burning-out 


Coats-Redfem’s models and activation energy is found to decrease 
with increasing palm shell composition in coal; however, a reverse 
trend is observed for the pre-exponential factor. 

While most of these studies about combustion of biomass and 
coal are conducted under a certain heating rate and usually 
selected biomasses mainly include wheat straw [13], corn straw 
and sawdust. Furthermore, kinetics analyses associated with the 
combustion are mainly focused on a certain stage and normally 
on individual biomass or coal [14-16], The use of bituminous coal 
with herbaceous biomass briquettes such as beetroot and switch- 
grass, and assessment about the effect of various parameters, like 
blending ratio and heating rate with wide scope, on co-combustion 
kinetics are still scarce. Little information is available to date on the 
kinetics model of co-combustion of biomass and coal based on 
two-stage scheme. This paper presents kinetics on co-combustion 
of bituminous coal and herbaceous biomass (beetroot and switch- 
grass) for a wide range of blending ratios and heating rates. The 
main purpose of the current study is to explore the effect of heat¬ 
ing rates, sample kinds, and biomass ratios on co-combustion 
kinetics based on two-stage scheme. Isoconversional and model¬ 
fitting methods for estimating kinetic parameters are compara¬ 
tively evaluated. The Arrhenius parameters are analyzed to assess 
if there are any compensating effects. 


2. Methods 

2.1. Experimental facility and test samples 

A thermogravimetric analyzer (TGA/SDTA851, Mettler Toledo) 
with a precision of 0.001 mg was used to conduct the co-combustion 
experiments. The beetroot and switchgrass residues were provided 
by University of Oakland, and Chinese pulverized bituminous coal 
were selected. The mixtures were obtained from a platinum crucible 
in which the beetroot/bituminous coal and switchgrass/bituminous 


coal blends were homogenized with biomass weight percentages of 
20,40,60 and 80 wt.% by using a 180 pm sieve to sieve the samples. 
Prior to thermogravimetric tests, biomass and coal were well-mixed 
mechanically and pre-dried for 10 h at 373 K. The kinetic character¬ 
istics of the overall process are determined by the velocity of 
advance of this interface into the unchanged reactant and the varia¬ 
tion in its effective area with time. Ideally, kinetic studies would be 
concerned with single crystals of simple and constant geometry but, 
in practice, it is often necessary to use assemblages of crystallites 
which usually contain a range of sizes. No unifying theoretical con¬ 
cepts have been recognized which can be used to provide satisfac¬ 
tory criteria for the comprehensive classification of the kinetics 
and mechanisms of reactions involving solids [17], The rate deter¬ 
mining step in any solid phase reaction can be either (i) diffusion, 
i.e. the transportation of participants to, or from, a zone of preferred 
reaction, or (ii) a chemical reaction, i.e. one or more bond redistribu¬ 
tion steps, generally occurring at a reaction interface. In mechanistic 
studies, measurements are frequently concerned with rate pro¬ 
cesses from which the effects of diffusion are specifically excluded. 
The initial mass of each sample was maintained at 10 mg ± 0.5 mg. 
The detailed method had been published from my previous work 
in the Refs. [18,19], The proximate analysis and ultimate analysis 
of the parent samples were listed in Table 1. The combustion pro¬ 
cesses were performed with heating rates of 10, 30, 50, 70 and 
90 °C/min, respectively. To ensure reproducibility, the experiments 
were repeated three times. The results indicate that a good 

Table 2 

Uncertainties of the measured parameters and calcu¬ 
lated parameters. 

Parameters Uncertainty (%) 

Temperature ±0.01 

Mass ±0.003 

Apparent activation energy ±0.01 
Pre-exponential factor ±0.014 


Table 1 

Proximate and ultimate analysis of samples. 


60B40C 

40B60C 

20B80C 


60S40C 

40S60C 

20S80C 



HHV (MJ/kg) 


17.01 

16.36 

16.72 

17.86 

18.50 
19.43 
16.45 
15.70 
15.29 
16.68 

16.50 
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reproducibility is maintained for each run because the relative 
deviation was generally within ±1.0%. 

The uncertainties of the measurement in the experiment are 
dependent on the experimental conditions and the measurement 
instruments. Therefore, the uncertainties of independent 


parameters are firstly computed, and consequently uncertainties 
of dependent parameters are calculated based on their relationship 
with independent parameters described in the Section 2.2. The 
uncertainties of the parameters in Table 2 were estimated by using 
the method in Refs. [20-22], 
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2.2. Kinetic analysis 


In order to calculate the kinetics much accurately, the main 
combustion region may be divided into two stages of pre-peak 
and post-peak by the maximum point in DTG. The reference value 
of activation energy of fuel around peak is calculated using the 
Ozawa-Flynn-Wall equation [23-28], which can be expressed as 
follows: 

,o s"= l “> ! [^)]- 2 ' 3,5 - 0457 & <’> 


where C(a ) is integral representation of the mechanism function, (I 
the heating rate (°C/min), T the absolute temperature (K), A the pre¬ 
exponential or frequency factor (s^ 1 ), E a the activation energy 
(J/mol), and R the universal gas constant [J/(mol K)], a the degree 
of conversion at time t and it can be calculated by the following 
equation 


where m is the weight of samples at any time, m 0 the initial weight 
of samples, and m x the weight of samples at the end of the burning 
process, respectively. By the way, the mass fraction of samples, lon¬ 
gitudinal coordinate in TG curves can be expressed as 


(3) 


By plotting log (j3) vs. 1/T, initial E a can be related to certain 
degree of conversion. 

The integral method based on the Coats and Redfern (CR) equa¬ 
tion [29-31] is a popular non-isothermal model-fitting method 
that requires an assumption be made regarding formula of the 
mechanism function g(a), which is expressed as follows: 


In [g(oc)/T 2 ] = In (AR/pE a ) - E a /RT 


(4) 


By substituting common solid reaction mechanism functions 
into Eq. (4) respectively [32], and plotting In [g(a)/r 2 J vs. 1/T, E a 
and A of fuel in each stage can be obtained. The mechanism func¬ 
tion will be selected solely by comparing the E a values with refer¬ 
ence value obtained from the Ozawa-Flynn-Wall equation. 


3. Results and discussion 


3.1. TG and DTG curves 


Although the TG and DTG curves of samples for a wide range of 
heating rates and blending ratios are obtained, for brevity, only 
some of the results are discussed. Fig. 1 illustrates the TG and 
DTG curves of 40B60C and 40S60C at different heating rates and 
all samples at heating rate of 10 °C/min. 

From Fig. 1, as heating rate increases, TG profiles shift to higher 
temperatures, and the rate of weight loss increases, which is indic¬ 
ative of the advance in combustion and complete combustion at 
lower temperatures. The mean combustion rate, and the maximum 
combustion rate of samples which mainly occurs at the volatile 
matter combustion region, increase with increasing the heating 
rate, as heightening of mass transfer results from the enhancement 
of heat transfer. As biomass blending ratio increases, TG profiles 
shift to lower temperatures, and the maximum combustion rate 
of samples which increases due to more volatile matter in biomass, 
while the combustion rate which occurs at fixed carbon combus¬ 
tion region decreases due to less fixed carbon in biomass. After 
the moisture release peak, pure beetroot and beetroot/coal blends 
showed 3-4 peaks for combustion of volatiles and fixed carbon, 
reflecting the existence of complex volatiles compositions for beet¬ 
root [18], The main combustion stages of beetroot, switchgrass and 


their blends with coal mainly focus at temperature range of 
200-400 °C due to beetroot and switchgrass being two typical 
herbaceous plants which are composed of more hemi-cellulose 
and cellulose. The biomass is mainly composed of hemi-cellulose, 
cellulose and lignin, which decompose at the temperature range 
of 225-325, 305-375 and 250-500 °C, respectively [33]. 


3.2. Kinetic parameters 

Fig. 1 reveals that there is a maximum value (DTG) at the tem¬ 
perature of 320 °C for 40% beetroot/60%bituminous coal blend at 
heating rate of 10°C/min. The linear fitting is shown in Fig. 2 by 
plotting log (£) vs. 1 IT for the blend (a = 0.33). R and SD are the cor¬ 
relation coefficient and the standard deviation of linear fitting, 
respectively [34,35], The reference value of activation energy of 
the blend is obtained, 143.50 kj/mol. The temperature range at 
pre-peak can be taken from 260 to 320 °C. Based on Eq. (4), by plot¬ 
ting In jg(a) /T 2 j vs. 1 IT as shown in Fig. 3, the dominant mechanism 
is fitted to the Avrami-Erofeev equation, G(a) = [- ln(l - a)] 4 , 
which is controlled by the random nucleation and then growing 
process. Nucleation is the process by which the reaction interface 
is initially established [17], at perhaps a limited number of points 
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tic parameters of beetroot/bitumir 




Samples 


20B80C 


40B60C 


60B40C 


80B20C 


Beetroot 


10 

30 

50 

70 

90 


420.00-480.00 1.64E+02 

465.00-515.00 3.45E+01 

496.00-546.00 1.17E+01 

496.00-546.00 4.25E+01 

524.00-574.00 7.00E+00 


53.63 

48.38 

43.80 

51.28 

41.54 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


0.01 

0.01 

0.01 

0.01 

0.01 


[l-(l-ar]/0.8 
[l _ (1 _ a) o-s] /0 .8 
[l-(l-a)°- 5 ]/0.5 
— [l — (1 — a )°' 8 ] /0-5 
-ln(l - a) 


10 

30 

50 

70 

90 


280.00-340.00 2.47E+05 123.96 


290.00-350.00 3.94E+03 

290.00-350.00 1.24E+04 

290.00-350.00 8.27E+03 

280.00-340.00 4.79E+04 


100.53 

105.97 

105.86 

105.92 


-0.99 


0.05 


-0.99 

-0.99 

-0.99 

-0.99 


0.03 

0.02 

0.01 

0.03 


[— In (1 - a)] 4 


10 

30 

50 

70 

90 


260.00-320.00 1.12E+06 

280.00-340.00 7.95E+05 

280.00-340.00 3.13E+06 

290.00-350.00 2.93E+06 

290.00-350.00 3.95E+07 


118.76 

114.85 

120.04 

117.50 

126.72 


10 

30 

50 

70 

90 


260.00-320.00 1.01E+08 142.53 

280.00-340.00 3.63E+08 149.27 

280.00-340.00 1.37E+08 156.73 

290.00-350.00 9.13E+07 158.97 

280.00-340.00 2.67E+08 162.61 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


0.01 

0.01 

0.02 

0.02 

0.04 


10 

30 

50 

70 

90 


260.00-320.00 4.18E+11 

240.00-300.00 8.75E+13 

260.00-320.00 7.73E+13 

260.00-320.00 1.65E+13 

260.00-320.00 2.96E+13 


162.31 

181.01 

185.12 

177.99 

181.34 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


[— ln(l - a)] 3 


10 


50 

70 

90 


260.00-320.00 2.37E+06 102.94 


280.00-340.00 3.17E+06 

290.00-350.00 8.07E+06 

290.00-350.00 1.10E+07 

300.00-360.00 2.55E+07 


102.27 

105.18 

105.86 

108.64 


-0.99 


0.01 


-0.99 

-0.99 

-0.99 

-0.99 


0.01 

0.02 

0.02 

0.02 


[-ln(l-a)] 3 


Table 4 

Combustion kinetic parameters of s 
Samples [1 (°C/mii 

20S80C 10 

30 
50 
70 
90 


ituminous coal at pre-peak. 

t(°C) _ Ms- 1 ) 

273.00-333.00 9.26E+12 

293.00-353.00 3.16E+13 

302.00-362.00 2.53E+13 

308.00-368.00 4.74E+15 

305.00-365.00 1.33E+12 


266.00-326.00 

295.00-355.00 

303.00-363.00 

304.00-364.00 

306.00-366.00 

272.00-332.00 

295.00-355.00 

303.00-363.00 

310.00-370.00 

313.00-373.00 


4.97E+12 

1.15E+12 

1.08E+12 

1.07E+13 

1.94E+08 

1.26E+08 

2.78E+08 

6.23E+08 


274.00-334.00 

287.00-347.00 

298.00-358.00 

307.00-367.00 

306.00-366.00 

269.00-329.00 

286.00-346.00 

297.00-357.00 


1.82E+09 

1.23E+10 

1.62E+10 

8.05E+09 

5.06E+09 

7.69E+15 


Eg (kj/mol) 


209.73 

213.08 

210.95 

239.41 

192.88 


181.31 

181.28 

190.99 

139.09 

135.51 
139.01 
142.18 
143.11 

143.52 


145.72 

154.36 

155.66 

152.60 

149.11 

211.24 

223.64 

163.70 

167.87 

160.68 


R 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.98 

-0.99 

-0.99 

-0.99 

-0.99 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


SD 


0.06 

0.04 

0.05 

0.07 

0.05 

0.06 

0.04 

0.05 

0.07 

0.05 

0.04 

0.03 

0.04 

0.10 

0.06 

0.06 

0.03 

0.04 

0.05 

0.04 

0.09 

0.05 

0.03 

0.04 

0.03 
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Samples 


60B40C 


80B20C 


Beetroot 


parameters of beetroot/bituminous coal at post-peak. 


P (°C/min) _tTO_ 

10 480.00-530.00 

30 515.00-575.00 

50 546.00-596.00 

70 546.00-606.00 

90 574.00-624.00 

10 340.00-400.00 

30 350.00-410.00 

50 350.00-410.00 

70 350.00-410.00 

90 340.00-400.00 

10 320.00-380.00 

30 340.00-400.00 

50 340.00-400.00 

70 350.00-410.00 

90 340.00-420.00 

10 320.00-380.00 

30 340.00-400.00 

50 340.00-420.00 

70 350.00-430.00 

90 340.00-420.00 

10 320.00-380.00 

30 300.00-360.00 

50 320.00-380.00 

70 320.00-380.00 

90 320.00-380.00 

10 320.00-400.00 

30 340.00-430.00 

50 350.00-430.00 

70 350.00-420.00 

90 360.00-440.00 


7.17E+07 
1.63E+02 
7.97E+01 
7.97E+01 
1.11E+01 
9.86E+04 
5.64E+02 
1.40E+02 
1.01E+02 
1.03E+03 
7.58E+03 
3.74E+02 
5.19E+02 
8.29E+01 
2.52E+01 
1.44E+04 
1.74E+03 
2.26E+01 
6.56E+00 
1.25E+01 
3.73E+04 
1.41E+08 
1.24E+08 
1.31E+08 
2.31 E+08 
3.88E+01 
2.03E+01 
3.00E+01 
1.44E+01 
4.78E+01 


Eg (kj/mol) 


127.36 

56.51 

54.20 

54.20 

43.96 


119.22 

90.65 
83.52 
80.93 
87.09 

94.55 

77.29 
77.37 
65.73 

55.66 
88.99 
75.82 

51.29 
43.64 

45.55 
89.27 
126.30 
125.08 
124.03 
127.15 


R 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.98 

-0.98 

-0.99 

-0.98 

-0.99 

-0.99 

-0.98 

-0.98 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.99 

-0.98 


SD 

0.01 

0.01 

0.01 

0.01 

0.01 

0.02 

0.02 

0.04 

0.04 

0.08 

0.07 

0.07 

0.09 

0.09 

0.07 

0.11 

0.11 

0.05 

0.03 

0.05 

0.13 

0.01 

0.06 

0.02 

0.01 

0.01 

0.02 

0.04 

0.03 

0.04 


-[l-(l-a)-°- 2 ]/0.2 

[l_ ( l_a)°- 2 ]/0.2 

-[l-(l-a)-°- 2 ]/0.2 

-[1 - (1 - a)- 02 ]/0.2 

-[l-(l-a)-° 2 ]/0.2 

[-ln(l-a)] 4 


[-ln(l-a)] 3 


[ ln(l - a)] 4 


vitchgrass/bitu 


: coal at post-peak. 


20S80C 


40S60C 


60S40C 


10 

30 

50 

70 

90 


333.00-393.00 2.20E+06 133.33 


353.00-413.00 5.78E+04 

362.00-432.00 4.02E+02 

368.00-438.00 3.76E+02 

365.00-445.00 6.66E+01 


-0.99 0.13 


-0.98 

-0.99 

-0.99 

-0.99 


0.13 

0.06 

0.08 

0.05 


10 

30 

50 

70 

90 


326.00-396.00 5.91 E+03 

355.00-425.00 4.04E+02 

363.00-433.00 1.45E+02 

364.00-444.00 2.80E+01 

366.00-456.00 6.41 E+00 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


0.09 

0.03 

0.06 

0.04 

0.02 


10 


50 

70 

90 


332.00-402.00 9.95E+02 

355.00-425.00 5.41 E+02 

363.00-433.00 4.83E+02 

370.00-450.00 2.13E+01 

373.00-453.00 2.37E+01 


-0.99 

-0.99 

-0.98 

-0.99 

-0.99 


0.05 

0.04 

0.07 

0.02 

0.04 


80S20C 


10 


50 

70 

90 


334.00-404.00 2.07E+03 

347.00-427.00 8.99E+03 

358.00-438.00 1.20E+03 

367.00-447.00 2.03E+02 

366.00-456.00 6.00E+01 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


0.05 

0.03 

0.04 

0.04 

0.03 


Switchgrass 10 

30 


50 

70 

90 


329.00-409.00 5.40E+03 

346.00-426.00 6.96E+04 

357.00-447.00 9.68E+03 

369.00-449.00 1.51E+03 

366.00-456.00 1.15E+03 


-0.99 

-0.99 

-0.99 

-0.99 

-0.99 


0.03 

0.03 

0.03 

0.03 

0.02 


in the reactant, usually at a surface and possibly at an imperfection. 
The initially generated germ nucleus, consisting of a few atoms only 
is generally regarded as unstable due to the high ratio of surface 


strain to volume. It may, however, develop into a growth nucleus, 
which increases in size through advance of the reactant-product 
interface into the bulk of the crystallite, so reducing the relative 
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(a) Pre-peak period 




E a (kJ-mol J ) 


(b) Post-peak period 


E a (kJ-mol *) 


Fig. 4. Kinetic compensation effect of co-combustion of biomass and bituminous coal. 


importance of the surface strain term. Therefore, the real activation 
energy and pre-exponential factor of the blend are gained, 
118.76 kj/mol and 1.12 x 10 s s _1 . The best correlations for two 
stages are obtained with the activation energy and reaction mech¬ 
anisms of all samples, as shown in Tables 3-6. The high correlation 
coefficient R values and low SD values indicate that the correspond¬ 
ing reaction model satisfactorily fits the experimental data. 

From the tables, the dominant mechanisms associated with 
combustion kinetics for beetroot and beetroot/coal blends at pre¬ 
peak and post-peak are described by the Avrami-Erofeev equations 
(n = 3 or n = 4), which are controlled by the random nuclear produc¬ 
ing and growing process. The mechanisms of combustion for 
switchgrass and switchgrass/coal blends at pre-peak are best fitted 
by the Avrami-Erofeev equation, Z-L-T equation or Anti-Jander 
equation (three dimensional diffusion); whereas their mechanisms 
are described by the Avrami-Erofeev equation (n = 4) at post-peak. 
Diffusion-controlled solid-state reaction kinetics in a sphere, where 
diffusion in all three directions is all-important. In a diffusion-con- 
trolled reaction, numerous chemical reactions or micro-structural 
changes in solids take place through solid-state diffusion, i.e. the 
movement and transport of gas molecules in the solid phase [36], 


However, the mechanisms of combustion for blends are much 
different that of bituminous coal, which chemical reactions with 
different orders (n = 0.2, 0.5, 1, 1.5 for pre-peak; n = 0.8, 1.2 for 
post-peak) are presented. The chemical reaction order model 
implies that the rate-determining step is the chemical reaction 
which can be described by simple reaction order. The activation 
energies of co-combustion for beetroot and coal during the pre-peak 
period and post-peak period range from 100.53 to 185.12 kj/mol, 
and from 43.64 to 127.15 kj/mol, respectively. While for 
switchgrass/coal blends, their activation energies at pre-peak and 
post-peak range from 135.51 to 239.41 kj/mol, and 45.86 to 
133.33 kj/mol, respectively. The relatively high activation energies 
may reflect to some extent dependence of pore structure on tem¬ 
perature. This implies that for the burning stage, higher tempera¬ 
ture is required to provide the forward burning reaction. In 
general, reactions with higher activation energies are more temper¬ 
ature-sensitive [3], Activation energies of combustion for poplar 
wood, and rice husk were calculated using three different isocon- 
versional methods, Flynn-Wall-Ozawa, Kissinger, and ASTM I-1I 
[28], Quang-Vu et al. revealed that the combustion activation 
energy of hemicellulose varied from 103.8 to 44.8 kj/mol for the 
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spruce wood, and from 144.7 to 41.3 kj/mol for the birch wood at a 
10 °C/min heating rate based on the differential method [35], That 
of char was from 183.1 to 109.4 kj/mol for the spruce and from 
222.0 to 132.3 kj/mol for the birch. The activation energy of the cel¬ 
lulose was from 221.5 to 239.0 kj/mol for the spruce, and from 
204.7 to 236.7 kj/mol for the birch. Oyedun et al. [37] obtained 
the pyrolysis activation energy of plastics and biomass blends at 
only a 10 °C/min heating rate based on the Coats-Redfern method. 
Combustion activation energy of microalgae in two stages at a 
10°C/min heating rate also based on the Coats-Redfern method 
was from 63 to 81 kj/mol and from 113 to 126 kj/mol, respectively 
[38], Gil et al. [6] evaluated the combustion kinetics of bituminous 
coal, pine sawdust and their blends at a 15 °C/min heating rate 
based on the Coats-Redfern method. In the above literatures, only 
a single heating rate and several mechanism functions were consid¬ 
ered. These models although based on experimental data appar¬ 
ently not impugned by heat and mass transfer effects, are not 
truly valid over widely variable conditions [13]. 

The kinetics characteristics of mixture samples had no linear 
relation with blend ratios due to interactions between biomass 
and coal. 

The biomass ratios and heating ratios have no effect on mecha¬ 
nisms of beetroot/coal blends at pre-peak and post-peak, and 
switchgrass/coal blends at post-peak, but they have certain effect 
on the mechanisms of switchgrass/coal blends at pre-peak. Accord¬ 
ing to the analysis principle of variance [39], the blending ratios 
have significant influence (P<0.05) on the apparent activation 
energy of co-combustion for beetroot/coal and switchgrass/coal 
blends at pre-peak and at post-peak, while the heating rates have 
more significant influence (P< 0.05, higher F value) on the appar¬ 
ent activation energy of co-combustion for switchgrass/coal blends 
at post-peak than blending ratios. 

3.3. Compensation effects 

Moreover, by analyzing the data in above tables, a linear rela¬ 
tion between E a and InA of each sample is revealed within heating 
rate of 10 °C/min to 90 °C/min. Fig. 4 summarizes the results. There 
are much less correlation coefficients for 60B40C blend at pre-peak 
and pure beetroot at post-peak, because their activation energies 
have not obvious change with different heating rates. 

Both the pre-exponential factor and the activation energy 
increases or decreases with the variation of heating rate, so the 
two parameters can compensate their effect on the co-combustion 
rate according to the Arrhenius principle [40], The increase of the 
pre-exponential factor seems to attenuate the decrease of reaction 
rate constant caused by increase of the apparent activation energy 
with the variation of heating rate. Compensation effect makes it 
more complicated for co-combustion of coal and biomass in differ¬ 
ent heating rates, because the decrease of the apparent activation 
energy and the pre-exponential factor play opposite roles on com¬ 
bustion rate. The kinetics compensation effects are also discovered 
in biomass pyrolysis, and pyrolysis and combustion characteristics 
of coal chars, respectively [41,42], A marked compensation effect is 
presented for the gasification reaction of biomass residues [43], 

4. Conclusions 

The beetroot and beetroot/coal blends showed 3-4 peaks for 
combustion of volatiles and fixed carbon. The blending ratios had 
significant effect on the apparent activation energy of co-combustion 
for beetroot/bituminous coal or switchgrass/bituminous coal. The 
heating rates have more significant influence on the apparent activa¬ 
tion energy of co-combustion for switchgrass/coal blends at post¬ 
peak. The activation energies of co-combustion for beetroot and coal 


during the pre-peak period and post-peak period ranged from 100.53 
to 185.12 kj/mol, and from 43.64 to 127.15 kj/mol, respectively. 
While for switchgrass/coal blends, their activation energies at pre¬ 
peak and post-peak ranged from 135.51 to 239.41 kj/mol, and 
45.86 to 133.33 kj/mol, respectively. A linear relationship between 
In A and E a provides a way to predict the kinetic parameters. 
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